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Abstract. In this paper we extend the recent work of C.A. Braumann [l] to the case of 
stochastic differential equation with random coefficients. Furthermore, the relationship of the 
Ito-Stratonovich stochastic calculus to studies of random population growth is also explained. 

1. Introduction 

As mentioned in the paper by Carlos A. Braumann [1], there is the issue of the Ito-Stratonovich 
controversy. That is, the issue of which stochastic calculus, Ito or Stratonovich, to adopt in 
the SDE models of the population growth. It has been shown that Ito and Stratonovich 
calculus give different results and do not yield the same solutions to the SDE models, leading 
to this controversy on which calculus is more appropriate when modelling population growth 
in a random environment. Hence, creating an obstacle on the use of this modelling approach. 
Carlos A. Braumann clears up the confusion concerning this controversy by showing that the 
apparent difference between the Ito and Stratonovich calculus is due to the confusion based 
on the assumption that both Ito and Stratonovich employ the same type of mean rates, i.e. 
interpreting the mean rate as an unspecified "average" per capita growth rate. In fact, Ito and 
Stratonovich calculus will yield exactly the same results when coupled with the appropriate mean 
rate. It is proven that, when using Ito calculus, b{N) is the arithmetic average growth rate 
Ra{x, t), and when using Stratonovich calculus, b{N) is the geometric average growth rate Rg{x, t). 

2. The Model 

Let N = N{t) denote the population size (number of individuals, density) at time t of a closed 
population (no migrations) and assume that the initial population size A^(0) — Nq > is known. 
In a randomly varying environment, we shall refer to ^ as the total population growth rate and 
to the per capita growth rate, jf^, simply by the growth rate. 

We can model the dynamics by assuming that the growth rate -i- ^ is the sum of an "average" 
growth rate b{N) and perturbations caused by random environmental fluctuations. We can 
approximate these perturbations by a white noise a-^{t). Where the growth rate is for 

some > 0, some density-dependent function b(N) (which is the growth rate of the population) 
having a continuous derivative. 

We can model the dynamics using the stochastic differential equation (SDE): 

^^=H^) + '^?(t), A(0) = iVo>0, (2.1) 

where we write the per capita growth rate as an "average" density-dependent rate 6(A) per- 

turbed by a white noise (J^{t), to take into account the effect of random environmental fluctuations. 

The b{N) represents the "average" growth rate in the population size. The (j£,{t) represents 
the random and uncertain movement in the population size which can be attributed to random 
environmental fluctuations that perturb the per capita growth rate. This perturbation is assumed 
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to be a stationary stochastic process and can be reasonably approximated by a white noise process. 
a(t) > is the volatihty and can be regarded as adding noise or variabiUty to the fluctuations in 
the population size, and £,{t) is a standard Gaussian white noise process. Furthermore, ^{t) is a 
generalised derivative of the standard Brownian motion process W{t) and is therefore equal to 

dWjt) 
dt ■ 

The SDE given in (2.1) can be rewritten as the following: 



dN 

— = b{t, u;)N{t) + a{t, iu)Nmt). (2.2) 

Here, lu Q represents the random environmental scenario (event) in the set D, of all possible 
environmental scenarios, on the probability space structure, (ri,jF, P). By a scenario lo, we mean 
a specific combination of environmental conditions that a population might be subjected to. 
We can rewrite this differential equation (2.2) in terms of a Brownian motion process, where we 
simply substitute £,{t)dt with dWt, since /p ^{s)ds — W{t) — /p dW{s). Hence, the basic model of 
the population growth in a random environment is given by: 

dN{t) ^ N{t)\b{t,Lu)dt + a{t,uj)dWt], N{0) ^ No > 0, (2.3) 

or in integral form 

N{t) = No+[ N{s)b{s)ds+ [ N{s)a{s)dWs, (2.4) 
Jo Jo 

where is a standard one-dimensional Brownian motion. b{-,Lu) and cr(-,a;) are assumed to be 
adapted and satisfy the integrability condition: 

T 

{\b{t)\ +a{tf)dt < oo 

almost surely, for every T £ (0, oo). This integrability condition ensures that N{t) < oo Vt <E M+. 

Since N{t) is continuous, the first integral in Equation (2.4), can be defined as a Riemann integral. 
However, problems now arise with the definition of the second integral of Equation (2.4), in 
which W{t) oscillates too rapidly to be defined in the usual Riemann-Stieljes sense (which follows 
ordinary calculus rules). The second integral in Equation (2.4) contains the Brownian motion 
component which defines this integral as a stochastic integral. This integral cannot usually be 
defined as a classical Riemann-Stieltjes integral due to the fact that the limits of Riemann-Stieltjes 
sums differ according to the choice of intermediate points in the integrand function. As a result, 
there are many alternative definitions of these stochastic integrals according to the choice of such 
intermediate points. 

The most commonly used integrals in the literature are the Ito and Stratonovich integrals. The 
Ito integral has nice probabilistic properties, which includes, apart from being a martingale, the 
property of zero expectation as well as having a convenient expression for its variance. However, 
it follows non-ordinary calculus rules. Stratonovich calculus, on the other hand, follows ordinary 
calculus rules. We will examine the problems that arise in the interpretation of (2.3) when the 
SDE is taken in the Ito and Stratonovich sense. 







In the next two sections, we will introduce the important concepts of the Ito and Stratonovich 
calculus respectively. Thereafter, in Section 5, we will discuss and represent the relationship 
between Ito and Stratonovich calculus. In Section 6, the controversy itself will be discussed along 
with its resolution. 
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3. Ito Calculus 
In Ito calculus we can express the Ito SDE in the following form: 

dXt = fi{Xt,t)dt + a{Xt,t)dWt, 
or equivalently in integral form as 



Xt^Xo+ f fi{Xs,s)ds+ f a{Xs,s)dWs 

JQ Jo 



Definition 3.1. The Ito Integral 

Suppose that W{t) is a Brownian motion process and that X{t) is a stochastic process. Consider 
a partition of [0,T], = to < ti < ■ ■ ■ < tn = T , then the ltd integral of X w.r.t. W is a random 
variable 



f XtdWt lim y2Xt^{W{t^+i) - W{t^)). 

""J j=0 



Notice, in the summation, the function X is defined at the left-hand point, i.e. the value of X at 
the beginning of each timestep is used, this is of crucial importance. 

Tiieorem 3.2. Ito's Lemma 

Let X{t) be a generalised Brownian motion process or an ltd process. That is, let X{t) have the 
following dynamics 

dXit) ^ a{Xt,t)dt + b{Xt,t)dWt, 
where Wt is a Brownian motion process. 

Let F{Xt,t) be a function with continuous second derivatives, where F and X have a functional 
dependence. 

Then F{Xt,t) is also an Ltd process and has the following dynamics 

dF BF 1 B'^F 

dF{Xt,t) = —{x,t)dt+—ix,t)dXit) + - — {x,t){dXit)f 

BF BF 1 B^F \ BF 

— {x,t)+a{x,t) — {x,t) + -b^{x,t) — {x,t)\ dt + b{x,t) — {x,t)dWt. 

Hence, F is also an Ito process, but with adjusted drift rate given by ^{x, t) + a{x, t)^{x, t) + 
^b'^(x,t)^^{x,t) and a scaled variance, b^{x,t)^{x,t). 

We may use Ito Lemma to solve (2.3), in the form 

d\iiN{t) ^-fit)dt + a{t)dWt, N{0) = No > 0, 

where 

7(i):=6(t)-ia2(i), (3.1) 

or equivalently: 

N{t) = Nn exp |^ j{s)ds + a{s)dW, | . (3.2) 
We note from (3.2) that N{t) > for alH > provide that iVo > 0. 

We shall refer to the quantity of (3.1) as the rate of growth of the population N, because of the 
a.s. relationship 

^lim^-(lnAr(T)-y^ 7(s)ds) = 0, (3.3) 

valid when the variance a(-) — cr'^(-) is bounded, uniformly in {t,uj); this follows from the strong 
law of large numbers and from the representation of (local) martingales as time-changed Brownian 
motions. 
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4. Stratonovich Calculus 

In probability theory, the Stratonovich integral is a stochastic integral, the most common al- 
ternative to the Ito integral. The appeal of Stratonovich calciilus is that in certain circumstances, 
integrals in the Stratonovich definition are easier to manipulate. Unlike the ltd integral counter- 
part, it is defined such that the chain rule of ordinary calculus holds for the stochastic integrals. 
Perhaps the most common situation in which these arc encoimtcrcd is as the solution to SDEs. 
These Stratonovich SDEs are equivalent to Ito SDEs, apart from the notation a{Xt,t) o dWt, 
where the " o " simply indicates that we are working in the Stratonovich sense. Furthermore, it 
is possible to convert between the two whenever one definition is more convenient for our purposes. 



In Stratonovich calculus we can express the Stratonovich SDE in the following form 

dXt = fJ,{Xt,t)dt + a{Xt,t) o dWt, 
or equivalently in integral form as 



Xt=Xo+ [ n{X„s)ds+ I a{X„s)odW,. 
Jo Jo 

Definition 4.1. The Stratonovich Integral 

It is defined in a similar manner to the Riemann integral, i. e. as a limit of Riemann sums. 
Suppose that W{t) is a Brownian motion process and that X{t) is a stochastic process. Consider 

a partition of [0,T], = to < ti < ■ ■ ■ < tn = T, then the Stratonovich integral of X w.r.t. W is 
a random variable 

pT n-l 

/ Xt o dWt := hm - W{ti)). 

Here, the function X is evaluated in the middle of each timestep (i.e. choose value of process at 
midpoint of each subinterval) . In the definition of the Ito integral, the same procedure is used 
except for choosing the value of the process X at the left-hand point of each subinterval, i.e. Xt^ 
in place of X(_t^_^^+t.y2- 

In a similar fashion to Ito's Lemma, Stratonovich satisfies 

/dF OF \ dF 

dF{Xt,t) = l—{x,t) + a{x,t) — {x,t)j dt + b{x,t) — {x,t)dWt. 



5. The Ito-Stratonovich Relationship 
Theorem 5.1. Conversion Formula 

Let X be a stochastic process, in particular, an R6 process satisfying the SDE, 
dXt = a{Xt,t)dt + a{Xt,t)dWt. Let a be some function of X and t. Furthermore, to 
make apparent the distinction between the two integrals, we will adopt the subscript I to indicate 
an ltd integral and the subscript S to indicate a Stratonovich integral. 



Then conversion between the R6 and Stratonovich integrals may be performed using the formula 



1 f'^ da 



as{Xut)dWt^ I ai{Xt,t)dWt + - —{x,t)a{Xt,t)dt. (5.1) 



Proof. 



f as{Xt,t)dWt- [ aiiXt,t)dWt = 
Jo Jo 

.^Jg ( ^fe)^^.) ,,) (£M±£M,.)) _ 
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By the Mean Value Theorem, we obtain 



Q 



3=0 



-X{t,) {Wit,+,)~Wit,)) 



= 1 1™ 'Y.^{x,t){x{t,+^) - x(tMW{t3+i) - Wit,)) 



2 n- 



3=0 
n-1 



dx 



3=0 
n-1 



da 

dx 



= 1 1™ y2^{X,t)AXtAWt 



3=0 

1 d 



3=0 

Which is the Ricniann sum of 

1 da, , , , 

Giving us the required result. 

From Theorem 5.1, it is evident that the Ito SDE, 

dX{t) ^ a{Xt,t)dt + a{Xt,t)dWt 
is equivalent to the Stratonovich SDE, 



dX{t) 

and that the Stratonovich SDE 



a{Xt,t)-^aiXt,t)^{Xt,t) 



dt + a{Xt,t)dWt 



dX(t) ^ a{Xt,t)dt + P{Xt,t)dWt 



is equivalent to the Ito SDE 
dX(t) = 



a{Xt,t) + ^P{Xt,t)^{Xt,t) 



dt + P{Xt,t)dWt. 



In the population dynamics context, the Stratonovich SDE, given by 



dN{t) = b{t)N{t)dt + a{t)N{t)dW{t). 

Using the Ito-Stratonovich conversion formula (5.5), we have a{Nt,t) = b{t)N{t) and 
a{t)N{t), where ^{Nt,t) = a{t). Hence, (5.6) is equivalent to the Ito SDE 



1 



b{t)N{t)--a{Nut)a{t) 
b{t)N{t) - \a\t)N{t) 



dN{t) = b{t)N{t) - ^aiNt,t)a{t) dt + a{Nu t)N{t)dW{t) 

dt + a{t)N(t)dWt. 
In a similar fashion we can obtain the reverse conversion formula. 
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6. The Controversy 

Many qualitative differences have been uncovered between Ito and Stratonovich calculus. In 

particular, there are instances in which Stratonovich calculus predicts, for the population, 
non-extinction and the existence of a stochastic equilibrium, whereas, at the same time, Ito 
calculus will predict population extinction. So, it seems, which calculus one uses does have 
important consequences. This fact has resulted in there being much controversy over which 
calculus is more appropriate to employ when finding a solution to the SDE. 

Considering the dramatic differences in predictions concerning important issues like extinction, 
which calculus should one trust? This is a major obstacle to the use of these stochastic models. 

Braumann (2003) resolved the issue of the Ito-Stratonovich controversy for the density- 
independent growth model, where b{N) = b is identically constant, in a random population 
environment. Braumann then extended these results in a random environment for the general 
density-dependent population growth model. It is revealed that the possible reason for this 
controversy, is the subtle fact that the same per capita "average" growth rate, '5' is used in 
both the Ito and Stratonovich calculus. Therefore, the issue here regards the meaning and 
interpretation of this average, since it is not elucidated what type of "average" is being referred 
to. Furthermore, it is inherently assumed that both the Ito and Stratonovich calculus make use 
of the same average, this is of course an incorrect assumption! 

Hence, this issue of the "average" needs to be addressed and clarification needs to be made of 
what type of average each method uses. In fact, it is found that the interpretation of b is different 
when considering population dynamics. When one decides to use Ito calculus in obtaining a 
solution to the SDE, b is interpreted as the arithmetic average growth rate. However, if ones 
chooses to implement Stratonovich calculus, b is interpreted as the geometric average growth 
rate. The differences between these two types of averages results in the dramatic differences 
between the Ito and Stratonovich calculus to disappear, yielding exactly the same solutions in 
both instances. 



Thus the differences are merely due to the absence of clarification of the meaning of b. So, all 
that is required is to match the appropriate average with the correct type of calculus, and exact 
same results will be obtained, putting to rest the Ito-Stratonovich controversy. 

6.1. Types of Averages. We will denote by ]Et,x['] := lE[-|A'^(t) = x], as the expectation condi- 
tioned on the knowledge that at time t the population size N{t) is x. 

6.1.1. The Arithmetic Average. The first type of average we have already mentioned is the arith- 
metic average. It is simply given by the usual expected value, conditioned on the knowledge that 

N{t) = X, and given by i = ^[X]. 

Definition 6.1. We define the arithmetic average growth rate at time t, when the population size 
at time t is x, as 

R.ix,t) ,^l,^^.^mt + m-^_ (6.1) 
^ ' X AtlO At ^ ' 

Theorem 6.2. Let Ra{x,t) be the arithmetic average growth rate as defined above, we have equiv- 
alently 

Ra{x,t)=Et,^[b{t)]. (6.2) 



Proof. We have 



Ut+At rt+At 1 

j{s)ds + a{s)dWs \ . 
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It follows 



1 E,,^[N{t + At)]-x 

1 Et,,[a;exp{/.*+^* j{s)ds + /.*+^* a(s)(it^J - x] 
= — lim 

X Atin At 

^ Et,,[(exp{/;+^* b{s)ds} - 1) X exp{/;+^* a{s)dW, - \ a\s)ds]] 

Atj.0 At 

□ 

6.1.2. The Geometric Average. 

Definition 6.3. The geometric mean of a positive random variable X is defined as e'^I'"^"^^! 
Since 



P'geom — | 



\i=l 



= gln[(n^.)-] 
^ gE[ln(X)] 



Definition 6.4. We define the geometric average growth rate at time t, when the population size 
at time t is x, as 

1 ^.^ exp(Et,^[lnjV(^ + At)]) - x 

X AtiO At 



R,{x,t) := ^ lim --^V-ML----n'-^--;j; (g_3) 



Hence, the geometric average is obtained by transforming the quantities to be averaged to log 
scale then taking an ordinary arithmetic average and then revert to the initial scale by inverting 
the algorithm. 



Proposition 6.5. Let Rg{x,t) be the geometric average growth rate as defined above, we have 
equivalently 

¥.t,xVnN{t + A t)] -\nx 
\t[o At 



R,{x,t) = -ML^^^-vv.-r^^.,j-x.... _ 



Proof. This follows from the fact that when z ^ we have (e^ — l)/z ^ 1. We apply this result 
to z := Et^x[ln{N{t + At)/x)] which tends to when At i as follows: 

li^ E,,4HNit + At)/x)] ^ ^ 

AtlO At AtlO At 

(e^_l) z 

= iim X — — 

Atio z At 

= lim(^ 

AtiO At 

1 ,. (xe^ — x) 
= — lim 

X AtiO At 

= Rg{x,t). 

□ 
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Theorem 6.6. Let Rg{x,t) be the geometric average growth rate as defined above, we have equiv- 
alently 

Rg{x,t) =Et.Al{t)]- (6.5) 
where "f(t) is as defined in Equation (S.l). 

Proof. We have 

{pt+At pt+At ^ 

l{s)ds + cj{s)dWs \ . 

By Proposition 6.5 above we have 

E,,.[lnjVft + A^)]-lnx 

= hm 

AtiO At 

ft+At 

hm Kt X 



□ 



Corollary 6.7. When -/{t,uj) = b{Nt) ~ . 

Let Rg{x) be the geometric average growth rate as defined above, we have equivalently 

Rg{x,t)=b{x)-^a^. 

It seems as though the arithmetic and geometric averages are equivalent, where the geometric 
average substitutes the process N{t) by the process \nN{t). These definitions give an indication 
on how these two rates can be estimated from observed data. For example, to determine an 
estimate of Ra{x, t), we look at all the instances t for which N{t) is close to x and then we take 
the average of those N{t + At) as an approximation of E,t.x[N{t + At)]. 

6.1.3. Other Types of Averages. These two averages are not the only averages we can consider, 
there are many other types of averages that are possible, so when we refer to an "average" it is of 
extreme importance to specify which particular average is being referred to. Some other possible 
averages are the: harmonic, median and quadratic averages. 



Median fi — med{xi, . . . , a;„) 
Quadratic A = (i ELi ^^V" - (E[^'])^ 



6.2. Density-Independent Growth. Recall that there are two main ways in which to interpret 
SDEs: Ito calculus and Stratonovich calculus. They usually lead to different solutions. Let us 
first consider the density- independent growth rate model in a random environment, this model 
corresponds to a constant "average" growth rate, b{N) = b, 

dN{t) = bNdt + a{t)N{t)dW{t). 

To make the distinction between the two approaches more apparent, we will use the notation bj, 
to denote the "average" growth rate under Ito calculus, and bs to denote the "average" growth 
rate under Stratonovich calculus. 
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6.2.1. ltd Model. Let us consider the density-independent Ito calculus model: 

dN{t) = biNdt + a{t)N{t)dW{t). (6.6) 

It is convenient to work in the logarithmic scale by making the change of variables, Y{t) = In N{t), 
y — Inx. 

Applying Ito's Lemma to Equation (6.6), where F{N{t)) ~ liiN{t), we obtain 

dl.Nit).'-^dNit)^iq^idNit)r (6.7) 

^ dN{t) + UTn^)idN{t)f (6.8) 



I ~ -cr^it)] dt + a{t)dW{t). (6.9) 



N{t) ' ' 2 \N{t)^ 
1 
2 

Therefore, Y{t) = IniV(i), satisfies the SDE in Equation (6.6). 

This can alternatively be expressed in the equivalent integral form as 



Y{t) =Yo + (^bi -^<T^^ ds + adWs (6.10) 
1 



= Yo+\^bi--a'jt + aWt. (6.11) 

From Equation (6.11), we can conclude that Y{t) ^ N{Yo + [bi — \(j'^)t, a^t). From this we obtain 
a solution to Equation (6.6), which is represented as 



N{t) = No exp 



]t + aWt 



2 



Hence, N{t) has a lognormal distribution with expected value the knowledge that E[e°'^'] = es'^^*, 
given by: 

¥.[N{t)] = A^oexp[6/t]. 

From Equation (6.11), one obtains the asymptotic result Y{t) ~ (6/ — icr^)t ast +oo. Therefore, 
as i ^ oo, N{t) ^ oo or N(t) according to whether the "average" growth rate bj is larger 

2 2 

than ^ or smaller than 

6.2.2. Stratonovich Model. Let us consider the density-independent Stratonovich calculus model 

dNit) = bsNdt + a{t)N{t)dW{t). (6.12) 

Since Stratonovich calculus obeys the ordinary calculus rules, we have by the ordinary chain rule 
of differentiation 

dY{t) ^ d\nN{t) ^ ^^^^dNit) (6.13) 

= bsdt + (7dW{t). (6.15) 
Therefore, Y{t) = \nN{t) satisfies the SDE given in Equation (6.15). 

This can alternatively be expressed in the equivalent integral form as 

Y{t)=Yo+ f bsds+ f adWs. (6.16) 



From which one immediately obtains the solution 

Y{t) = Yo + bst + aWt. (6.17) 
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Since the integrand is constant, in this density- independent case, the Ito and Stratonovich integrals 
coincide. For both approaches, we have crWt — a{Wt — Wq) — aWt since Wq = 0. 
Since W{t) ^ N{0,t) (i.e. normally distributed with mean zero and variance t), we conclude that 
Y{t) ^ N{Yo + bs, a^t). From this we obtain a solution to Equation (6.12), which is given by 

N{t) = No exp [bst + aWt] . 
Hence, N{t) has a lognormal distribution with expected value given by: 



i[N{t)] = No exp 



rcr^ I t 



From Equation (6.17), one obtains the asymptotic result Y{t) 
t — > oo). Therefore, as t - 
according to whether the 



bst as t - 

■ oo, N{t) — > oo (i.e. growth without bound) or N{t) 
average" growth rate bs is positive or negative. 



oo (since ^ as 
(extinction) 



6.2.3. Conclusion. The long-term behaviour of N{t) for both interpretations of the SDE can be 
further analysed by examining the trajectory of N{t) in probability. Since, ^ — > a.s. when 
t — > oo, one easily notices that under Stratonovich calculus, N{t) oo when 6s > (probability 
of extinction is zero and there is a stochastic equilibrium) and N{t) when bs < (i.e. 
population extinction occurs with probability one). 

In a similar fashion, under Ito calculus, N{t) oo when bi > ^ and N{t) when 6/ < 

The differences between the Ito and Stratonovich approaches are now apparent. The behaviour 
appears to be different from the Stratonovich calculus. Hence, if one employs Ito instead of 
Stratonovich, the conditions for non-extinction and existence of a stochastic equilibrium are 
qualitatively different. This illustrates the consequences of the two approaches for the population 
behaviour. Using Stratonovich calculus, extinction would occur a.s. if the "average" growth rate 
6 < 0, but with Ito calculus, one can have extinction a.s. even for positive values of the "average" 

2 

growth rate 6 if 6 < ^ . 

Furthermore, the Ito calculus obtains different results compared to the deterministic model, this 
makes Ito calculus quite popular in modelling, hence, avoiding the issue of ignoring random 
environmental fluctuations. 



The approach taken here works for all density-dependent models and completely and exactly 
elucidates the difference between the two interpretations. It also exactly solves the problem of 
which calculus to use and how to use it. 



Ito SDE: N{t) = biNdt + aNdWt 

d\nN{t) = (^bi - iff^^ dt + adWt 

Stratonovich SDE: dN{t) ^ bsNdt + aNdWt 

d\nN{t) = bsdt + adWt 

6.2.4. Resolution oj which Average to use. 6.2.4.1 Ito 

Let us compute the two averages for the Ito SDE (6.6), we obtain from (6.17), 
Y{t + At) = lnx+ (bi - ^aA At + (j{W{t + At) - W{t)). 
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Therefore, Y{t + At) is normally distributed with mean Ina; + (6/ — ^)At and variance a^At. 
The conditional expectation is 

Et,^[Y{t + At)] = Et,SnN{t + At)] 
= lna;+ (^^^ ^ y 1 

Replacing into equation (6.3), we obtain 



1 exp(lnx+(bi-^)AO-x 
Ra{x,t) — — lim 

' ' X AtlQ At 



exp[(6/ - ^)At] - 1 

= lim 

AtJ,0 At 

^bi - ifT^ = bs- 

We also notice that N{t + At) = exjp{Y(t + At)) is lognormal with parameters Inx + (6/ — ^)At 
and a^At, and so its conditional expectation is 



Et,4N{t + At)] = Et^,[exp(y(t + At)] 



^2 



exp ( Inx + ( 6/ - —j At + a{W{t + At) - W{t)) 
exp [\nx+ [ hi - ^ \ At + ^cr^At 



= X exp (6/ At). 
Replacing into equation (6.1), we obtain 



„ / N 1 xexp(&/At)— X 
Ra(x,t) = - lim '- = bi. 

^ ' ^ X Atia At 



The conclusion is that 

Rg{x,t) = bi 2 , 
Ra[x,t) = bi. 

Hence, when using Ito calculus, the "average" growth rate bi is specified as the arithmetic 
average growth rate. 

6.2.4.2 Stratonovich 

Let us compute these two averages for the Stratonovich SDE model (6.12). Since N{t) = x, we 
obtain from (6.17), 

Y{t + At) = Ina; + bgAt + (T{W{t + At) ~ W{t)). 
Therefore, Y{t + At) ~ N{\n x + bsAt, a'^At) and so its conditional expectation is 

Et4Y{t + At)] = Et4lnN{t + At)] = Inx + bgAt. 
Replacing into equation (6.3), we obtain 
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1 cxp{lnx + bs At) - X 
R„{x,t) = — lim 

1 xe^^^' - X 

= — lim 

X Atio M 

1 , x{e^^^^ - 1) 

= — lim 

X Atio At 

= hs. 

N{t + At) = exp(F(t + At)) is lognormally distributed with parameters Ina; + bsAt and cr^At, 
hence its conditional expectation is 

Et,,[iV(t + At)] = Et,,[exp(y(t + At)] 

= Et^^[exp(ln2; + &sAt + a{W{t + At) - W{t)))] 

= exp ^Ina; + fe^At + ^cr^At 

= xexp [{bs + ]^(T^^ At 

where Et,^[e'^(^(*+^*)-^(*))] = e^-^'^*, since W(t + At) - W{t) - iV(0, At). 
Replacing this into equation (6.1), we obtain 

1 xexp((6s + ^a^)At)-x 
Ra(x,t) — — lim 

^ ^ X AUO At 

exp[(6s + ia2)At] - 1 
= lim f 

AtiO At 

= bs + ^f^^At. 

The conclusion is that 



Rg{x,t) = bs, 

Ra{x,t) = bs+^(7^. 

Hence, when using Stratonovich calculus, the "average" growth rate bs is specified as the 
geometric average growth rate. 

6.2.4.3 Conclusion 

The conclusion is that when using Ito calculus, the "average" growth rate, bj is specified as the 
arithmetic average growth rate, and when using Stratonovich calculus, the "average" growth 
rate, bs is specified as the geometric average growth rate 

This fact instructs us to replace the unspecified growth rate b by the specified average it truly 
represents. It is only in this manner that the results acquire meaning. Since bs is indeed the 
geometric average growth rate Rg, we can conclude that the solution of the Stratonovich SDE 
density-independent growth model is 



N{t) = No exp[Rgt + aW{t)]. 

Since 6/ is indeed the arithmetic average growth rate Ra, we can conclude that the solution of the 
Ito SDE density-independent growth model is 
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N{t) = No exp 



Ra 



— ]t + aW{t) 



or, alternatively it was shown that Ra — ^ = Rg, yielding 

N(t) ^ NQeyii>[Rgt + aW{t)]. (6.18) 

Hence, we can conclude that the two interpretations yield exactly the same solutions in terms 
of a specific average growth rate. Therefore, it does not matter which average we choose as 
long as it is clearly specified. With regard to the conditions under which extinction occurs, we 
conclude that both approaches predict population extinction or a stochastic equilibrium according 
to whether the geometric average growth rate is negative or positive. So, we can use either 
calculus indifferently as long as we are careful to use b for the appropriate average for that calculus. 



Once Ra{x,t) or Rg{x,t) — Ra{x,t) — ^ ' have been estimated, one can choose to es- 
timate bj{t) — Ra{x,t) — Rg{x,t) + ^ ' and use Ito calculus or choose to estimate 
bs{t) = Rg{x,t) = Ra{x,t) ^ and use Stratonovich calculus. 



It does not matter what choice one makes, the solution one obtains is the same. Therefore, the 
Ito and Stratonovich SDEs can be written in terms of these estimated quantities, Ra{x,t) and 
Rg{x,t): 

Ito SDE: - Raix, t) + a{t)£,{t) = Rg{x, + y + <mt) 



Stratonovich SDE: 
which are equivalent. 



6.3. Density-Dependent Growth. We can interpret (2.3) as an Ito SDE and it can be written 
as: 



dN{t)^N{t) bi{t)dt + a(t)dWt . (6.19) 
By the Ito-Stratonovich conversion formula given in Equation (5.3), where a{Nt,t) = b{t)N{t), 

da f 
dx ^ 



a{Nut) = a{t)N{t) and f?(iVt,t) = a{t). 



Therefore the Ito SDE is equivalent to the Stratonovich SDE 



dN{t) 



bi{t)N{t) - -(7'^{t)N{t) 



dt + a{t)N{t)dWt, 



dN(t) = N{t) 



bi{t) - 



aHt) 



dt + a{t)N{t)dWt 



= bsN{t)dt + a{t)N{t)dWt. 
This can be written in terms of the per capita average growth rate as 



1 dN{t) 
N dt 



bj{t) - 



a^it) 



a{t)m 



(6.20) 

(6.21) 
(6.22) 



2 

^bs + a{t)£,{t). 

This is similar to Equation (2.3) but simply interpreted as a Stratonovich SDE, where b{t) = bi{t) 

is replaced by b{t) ^ = bs- Whether we interpret equation (2.3) as an Ito or a Stratonovich 

SDE, the solution is a homogeneous diffusion process with diffusion coefficient (variance rate), 
N'^a'^ (which is the same in both the Ito and Stratonovich SDEs). The drift coefficient is, however, 
different; it is respectively for Ito and Stratonovich: 
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Ito 

Stratonovich 



^iiN) ^ N{t)bi{t) 



^l{N) = N{t) 



bi{t) 



a' it) 



To conclude, we have: 



Ito SDE 



Stratonovich SDE 



1 dN{t) 
N{t) dt 
dN{t) 
1 dN{t) 
N{t) dt 
dN{t) 



N{t)bsit). 

bi{t)+<j{t)m 

N(t)bi{t)dt + a{t)N{t)dW{t), 

bs{t) + <7{t)at) 

N(t)bsit)dt + a{t)N{t)dW{t). 



6.3.1. Deterministic Model. We need to clarify what the growth rate and "average" growth rate 
mean in terms of the observed population dynamics N{t). Let us consider the deterministic model, 
where a — 0, 



dN{t) = bd{t)N{t)dt. 



(6.23) 



We can define the growth rate (i.e. the per capita growth rate) at time t when the observed 
population size at time t is a;, i.e. N(t) = x, as 



bd{x) : 



1 dN{t) 



dt 



X 

— lim 

X AtiO 

i hm ^ 

X AtiO 



N{t + At) - N{t) 
At 

At) - a; 



The limit part represents the total growth rate, 



At 

dN{t) 
dt 



(6.24) 
(6.25) 
(6.26) 



at time t, which is then divided by 
the population size x to obtain the per capita growth rate. In this density-independent case, 
bd{x) does not depend on the population size x, however, we adopt the notation bd{x) since the 
dependence will be present in the more general density-dependence case. 

Alternatively, we could have obtained the solution of (6.23) as N{t) = Nocxp{bdt), where this 
gives N{t + At) — N{t) exp(6dAt), and substituting into equation (G.26) to obtain 

1 , xie''"^* - 1) 
— hm . 

X At^O At 



Furthermore, we know that as At | 0. 



At 



— tends to bd, to arrive at the same solution. 



However, in the stochastic model (2.3), i.e. cr ^ 0, N{t + At) is a random variable and it is then 
necessary to take some kind of average of iV(t-|- At), to obtain a possible estimate of the population 
size at time t -I- At. One possible way is to take the limit and determine the average afterwards, 
however, this is shown not to work since the limit itself is a generalised stochastic process and 
does not exist in the ordinary sense. So, instead, we follow the other approach by first taking 
the average and then computing the limit afterwards. But we have to be precise on what type of 
average we are using. 
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6.3.2. ltd Model. Let us see what these two averages turn out to be, under the ltd model. It is 
sometimes more convenient to work with the log transformed stochastic process Y[t) = In N(t) 
and with y = ln(a;). Using Ito's rule of calculus, one obtains, using the fact that in the limit as dt 
tends to {dt 0), (dt)^ = 0, {dt){dWt) = and {dWtf = dt: 

dN{t) = bi{t)N{t)dt + a{t)N{t)dW{t). (6.27) 

This Ito model has drift coefficient bi{t)N{t) and diffusion coefficient a^{t)N'^{t), where the growth 
rate bi(t) is given by; 

bi(t) = — lim — — i — L—t ^- = R (x.t). 

^ ^ X Atio At ^ ' ' 

The arithmetic average growth rate is then given by Ra{x,t) = bi(t). 

To compute the geometric average growth rate, let us consider the log scale Y{t) — \nN{t), then 
using Ito's lemma, we obtain the SDE (6.30): 

= ^ [Nit)biit)dt + a{t)N{t)dW{t)] + \ [a\t)N\t)dt\ (6.29) 

dt + (T{t)dW{t). (6.30) 



bi(t)-\a^{t) 



In terms of Y{t) = \nN{t), y = h\{x), the solution of equation (6.30) has drift coefficient bi{t) 
icr^(t) and diffusion coefficient cr^{t), where the growth rate bj{t) — ^a'^{t) is given by: 

bi{t) ~ -a it) ^ hm = R,{x, t). 

Therefore, the geometric average growth rate is given by 



Rg{x,t) = i^^^ — -^r^ ^ - Y^^^ ^ 

Hence, for the Ito SDE, we have that 

Ra{x,t) ^bi{t), 

2 

Rg{x,t)^blit)^^=bsit). 

These are respectively the arithmetic average growth rate (the expected value of average w.r.t the 
process N{t) as per the definition of the average (6.1)) and the geometric average growth rate 
defined by (6.3) for the solution of the Ito SDE. 

By the definition (6.3) of the geometric average, Rg{x,t) is the average w.r.t the process \nN{t), 
so by determining the dynamics of In N{t) we can obtain the drift rate which gives the per capita 
geometric growth average. 

To conclude bi{x,t) is the arithmetic average growth rate 

p . ,^ 1 KtAN{t + At)]-x 

Ra(x,t) := - hm — ; 

X Auo At 

and bi{t) ^ = bs is the geometric average growth rate 

Et,x[ In N{t + At)] -\nx 



Rg{x, t) := lim 



AtiO At 



Therefore, Rg{x,t) = Ra{x,t) ~ 



2 
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6.3.3. Stratonovich Model. For the Stratonovich model, we make use of an easier approach which 
is to convert the Stratonovich SDE to an equivalent Ito SDE using the Ito-Stratonovich conversion 
formula (5.5): 

Stratonovich SDE : ^^^^^ = bs{t) + cr(i)^(i), 

dN{t) = N{t)bsit)dt + <7{t)N{t)dW{t) 

equivalent Ito SDE : '^^J^^ = bsit) + + C7{t)^{t) 

dN{t) = N{t)bs{t)dt + N{t) {^-^^ dt + a{t)N{t)dW{t), 



Nit) 



dt + ait)Nit)dWit). 



2 

The solution of the Stratonovich SDE is a diffusion process with drift coefficient 

(bsit) + ^"2^^ ^(t) 8.nd diffusion coefficient a'^ (t) N'^ (t) , which is identical to the diffusion 
coefficient of the Ito SDE. 

If we now consider the transformation Y(t) — In Nit), which is instrumental in these deductions. 
Using Ito's Lemma we obtain 

= ,.ru\ dNit) + -^-lidNit))^, (6.31) 



dNit) ' ' 2 dNit) 



1 



Nit) 



Nit) { bsit) + ^ ) dt + ait)Nit)dWit) 



(6.32) 



= hsit)dt + ait)dWit). (6.33) 
The solution of equation (6.33) has drift bsit) and diffusion coefficient cr^(t). 

Alternatively, we could compute the geometric average growth rate by starting from the 
Stratonovich SDE and using the ordinary chain rule of differentiation, we obtain the SDE: 

^ [Nit)hsit)dt + (Tit)Nit)dWit)] 



Nit) 

= bsit)dt + ait)dWit). 

This can indifferently be interpreted as an Ito or a Stratonovich SDE. Since, the stochastic term 
has a constant coefficient, the correction term in the conversion method is now zero and the two 
approaches coincide. 

Therefore, using a similar reasoning as explained with the Ito SDE, for the Stratonovich SDE, we 
have 

Rgix,t)=bsit), 

= biit)-la' 

2 

Raix,t)^bsit) + Y 
= blit). 

These are respectively the geometric average growth rate and the arithmetic average growth rate 
for the solution of the Stratonovich SDE. This proves that the average used in Stratonovich calculus 
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is the geometric average and the average used in Ito calculus is the arithmetic average. Hence, 
from these results, we reach the following final conclusions: 
To conclude the arithmetic average growth rate is given by 



Ra[x, t) := - hm — = bsit) 

X AtiO At 



2 



Rg[x,t) := hm — bs[t). 



and the geometric average growth rate is given by 

Ef,^[ lniV(t + At)] - In a: 
ATIb At 

Again, for the Stratonovich SDE, we also have that Rg{x,t) — Ra{x,t) — 

7. Conclusion 

Under Ito calculus, we interpret the growth rate as the arithmetic average growth rate, 
Ra{x,t) defined by equation (6.f). 

Under Stratonovich calculus, we interpret the growth rate as the geometric average growth 
rate, Rg{x,t) defined by equation (6.3). 

Therefore, for Ito calculus, b{t) really means the arithmetic average growth rate, Ra{x,t) or 

2 

equivalently Rg{x,t) + ^ and for Stratonovich calculus, 6(t) really means the geometric average 
growth rate, Rg{x, t). 

It is shown, and finally concluded, in Braumann (2007) that both calculus lead to the exact same 
conclusions in terms of the conditions under which population extinction or the existence of a 

2 

stochastic equilibrium occur. Since if one takes into account the difference ^ between the two 
averages, the solutions under both approaches coincide. 

Hence, after clearing the confusion, Ito and Stratonovich calculus yield the same results. It is now 
possible to easily tackle the major obstacle to the use of these SDE models, however, care must 
be taken in using the appropriate type of average for each calculus If that care is indeed taken, 
Ito and Stratonovich will give the same results and draw the same conclusions. 
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